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Исследуются возможности решения задачи идентификации таких трудноопределимых физических парамет- 
ров математической модели участка магистрального газопровода, как удельный усреднённый по длине тру- 
бопровода коэффициент трения и фактическая средняя по трубопроводу температура при использовании 
изотермического варианта модели. Разрабатывается математический аппарат для такой идентификации на 
основе спланированного эксперимента. 
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Введение. Для решения задач автоматизации технологических процессов дальней транспорти- 
ровки природного газа характерна ситуация, когда используемая для управления или диагностики 
математическая модель процесса течения газа не является до конца параметрически определён- 
ной. Это касается ряда трудноопределимых констант, связанных с удельным сопротивлением тру- 
бопроводов, осреднением скорости потока по сечению, теплообменом их с окружающей средой, 
осреднением температуры по длине трубопровода при использовании изотермической модели и 
др. Поэтому аналитически построенную и, следовательно, структурно адекватную процессу мате- 
матическую модель необходимо идентифицировать параметрически. Это особенно важно для 
описания динамики магистральных газопроводов (МГП), т. к. посредством исследований установ- 
лено, что структура связи переменных состояния этого распределённого объекта такова, что его 
модель очень чувствительна к некоторым параметрам и имеет физический смысл только в опре- 
делённых и довольно узких диапазонах [1—3]. Исследования по определению этих диапазонов 
позволяют уточнить (параметрически идентифицировать) ряд трудноопределимых параметров, а 
также, опираясь на найденные диапазоны их возможных значений, решить задачу робастного вы- 
бора конструктивных и технологических параметров технологических участков МГП [3—6]. 

Формулировка задачи исследования. Для решения комплекса обозначенных здесь задач не- 
обходимо исследовать те аналитические возможности, которые предоставляют получаемые на 
стадии математического моделирования динамические характеристики распределённого объекта. 
При этом задачу целесообразно привязать к конкретному технологическому объекту — участку 
магистрального газопровода. Кроме того, её достаточно локализовать на оценке установившихся 
значений переменных процесса транспортировки газа, а также на отдельных конструктивных па- 
раметрах трубопровода. Необходимо привлечь и исследовать те взаимозависимости переменных 
состояния и параметров объекта, которые позволяет выделить математическая модель его уста- 
новившегося состояния, т. е. модель стационарного процесса течения газа по участку МГП [3—6]. 
Исходная математическая модель задачи исследования. Решение сформулированной вы- 
ше задачи опирается на математическую модель, характеризующую изотермический процесс! пе- 





1 Режим течения газа при интенсивном теплообмене с окружающей средой; принимается за парадигму для трубопроводов 
дальней транспортировки газа при его грунтовой прокладке. 
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редачи потоком газа массы вещества и преобразования его механической энергии, описывается 
математической системой следующего вида: 
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Модель (1) построена в базисе (р,9), где р(х,ё) — давление в произвольном сечении 


потока в произвольный момент времени (Па); 4 (х,Ё) — удельный (на единицу площади сечения 
трубопровода) массовый расход газа (кг / с.м?); В — универсальная газовая постоянная 


(Дж / моль .К ); Т — абсолютная его температура (К); и — молярная масса газа (кг / моль) ; 


р — диаметр трубопровода; ^, — удельный (на единицу площади боковой поверхности трубо- 
провода) коэффициент трения газа. Константы модели заданы следующими значениями: 
р = 1,22 м; К = 8,3144; р = 0,01604 моль, а ^, и Т подлежат идентификации. 

В связи с поставленной выше задачей решения проблем идентификации в рамках стацио- 
нарного состояния объекта исследования модель (1) необходимо представить структурно в не- 
сколько другом виде, обусловленном допущением о затухании всех переходных процессов. Оно 
связано с равенством нулю производных по времени всех переменных моделируемого процесса. 
Таким образом, необходима модель стационарного течения газа по трубопроводу. В базисе моде- 
ли (1) — (р,9) она имеет следующий вид: 


ф(х)_ 4.)..К.Т.р(х)-9(х) 
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Интегрирование уравнения (2 Ь) выявляет свойство инвариантности массового расхода к 
пространственной координате, которое моделируется условием Ух е [0,4 |9 (х) = 4% = с0П5ё . 


Тогда, интегрируя (2 а) при 4 (х) = 9, ‚ можно получить систему алгебраических уравнений, экви- 


валентную системе дифференциальных уравнений (2) 
9(х) =9; (а) 


2 рен [| 4х0 (Ь) я 
2.К.-Т- 49 В р 
Математическая модель решения задачи идентификации. Структура зависимости (3) по- 
зволяет получить явные выражения для температуры потока газа и коэффициента трения потока 
о стенки трубопровода. При необходимости их также можно преобразовать, для того чтобы в рас- 
чётах основываться не на фиксированном базисе, а использовать в качестве базисных перемен- 
ных те переменные состояния потока, значения которых известны и которые удобны для иссле- 
дователя. 

Для демонстрации возможности такого преобразования и корректности расчёта, основан- 
ного на нём, приведём вариант преобразования базиса относительно скорости потока газа. В этом 
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случае выражение (ЗЬ) примет вид (4) и будет содержать только одну переменную состояния по- 


тока — скорость: 
в.Т и? ЗЫ, 
1-9 + .х=0. 4 
2.и.и? | ". |= р (9) 


Из (4) можно выразить 7, (5) и Т (6): 
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Объединим (5) и (6) в систему. Чтобы избежать линейной зависимости, используем два 
набора опытных данных при различных начальных условиях (у, и, ии, и,). 
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Полученная ММ УМГП, разрешённая относительно температуры потока газа и коэффици- 
ента трения потока о стенки трубопровода, позволяет произвести идентификацию коэффициента 
трения потока газа о стенки трубопровода и средней температуры потока в ситуации, когда зна- 
чения этих констант либо не замеряются на всём протяжении участка, либо количество замеров 
не позволяет корректно рассчитать их усреднённые оценки. 

Хотя система (7) легко разрешима относительно неизвестных ). и Т , возникает проблема 
точности результата, поскольку погрешности в измерениях исходных величин могут привести к 
неточности итогового результата. Чтобы повысить точность идентификации, целесообразно ис- 


пользовать не два набора опытных данных, а п, т. е. и,, ии, при / =1,п . Построим всевозмож- 
ные сочетания этих измерений из п по 2 без повторений, разрешим для каждого из сочетаний 
систему (7) и вычислим усреднённую оценку /, и Т. С помощью этого, в соответствии с законом 
больших чисел, удастся получить искомый результат с высокой достоверностью. 


Построение комбинаций опытных данных для решения системы (7) позволяет увеличить 
выборку без увеличения числа опытов. В табл. 1 указаны соответствия числа измерений (п) чис- 


лу возможных сочетаний по два (С? ). 











Таблица 1 
Соответствие числа измерений числу сочетаний 
п 2 3 4 5 6 7 8 9 10 
С? 1 3 6 10 15 21 28 36 45 
































Пример решения задачи идентификации. Рассмотрим пример расчёта, произведённого со- 
гласно описанному методу. Исходные данные представлены в табл. 2. 

Простейшим вариантом будет являться расчёт на основе лишь двух наборов опытных дан- 
ных, необходимых и достаточных для однократного решения системы (7) и нахождения весьма 
приближённой оценки исследуемых констант. В более сложном варианте, т. е. при увеличении 
количества измерений, точность результатов возрастает. Поскольку этот случай представляет 
больший интерес, будем производить расчёты на основе 4 наборов опытных данных (измерений 
начальной и конечной скорости; табл. 3). 
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Таблица 2 
Исходные данные для идентификации 
Е, м р, м и, — 
МОЛЬ 
137000 1,195 0,016 
Таблица 3 
Начальная и конечная скорости 
№ и", м у, м 
1 13 81,0026 
2 12 29,2682 
3 11 20,066 
4 10 15,3994 




















В рассматриваемом случае мы ограничились 4 измерениями. Как видно из табл. 1, эти из- 


мерения могут быть использованы в 6 сочетаниях. 

Для каждой из комбинаций заданных измерений решим систему (7), найдём 6 пар значе- 
ний температуры и коэффициента трения. Путём вычисления среднего арифметического соответ- 
ствующих элементов этих пар найдём усреднённые оценки искомых величин. Результаты сведём в 


табл. 4. 


























Таблица 4 
Расчёт комбинаций и средних значений 
№ Изм. 1 Изм. 2 А, Т, м 
1 1 2 0,00096086 292,96 
2 1 Е 0,00096148 293,15 
3 1 4 0,00096131 293,091 
4 2 Е 0,00096349 293,76 
5 2 4 0,00096221 293,37 
6 Е 4 0,00096003 292,701 
Средн. — — 0,00096156 293,17 




















Результаты исследований. Таким образом, коэффициент трения идентифицируется значением 
0,00096156, температура — значением 293,17. Заметим, что реальное значение коэффициента 
трения при заданных параметрах 0,000961, а температуры — 293 К. Таким образом, точность 
идентификации констант при 4 наборах опытных данных оказалась достаточно высокой. 


























Таблица 5 
Квадраты отклонения от среднего значения 

№ (.-^) (т ег. 

1 0,491717 : 10-2 0,0459939 
2 0,00687775 : 10-2 0,000633233 
3 0,066689 10-2 0,00621534 
4 3,70436 ' 10-2 0,3454766 
5 0,421776 ' 10-2 0,03927496 
6 2,34615 ' 10-2 0,2185451 











Исследуем полученные результаты с точки зрения математической статистики, чтобы до- 
казать, что использование комбинаций позволяет на основе относительно небольшого числа из- 
мерений получить статистически достоверный результат. Вычислим квадраты отклонений Л иТ 


от средних значений для каждой из комбинаций (табл. 5). 
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Рассчитаем среднеквадратическое отклонение с, дисперсию ос”, а также стандартное от- 
клонение $ (оценку среднеквадратичного отклонения относительно математического ожидания 
на основе несмещённой оценки дисперсии) для ^. и Т (табл. 6). 


























Таблица 6 
Статистические и вероятностные показатели 
Показатель А т 
6 —\2 
У (^, -^) 7,037567 : 10-2 0,656139 
ТЕР 
5 1,083018 * 10-8 0,330691 
о? 1,172928 * 10-2 0,109357 
5 1,407513 ' 10-2 0,362254 
Среднее значение 0,00096156 293,17 
Относительная ошибка 0,0011263 0,0011280 











Заключение. Таким образом, теоретические расчёты и экспериментальная проверка показали, 
что предложенный метод позволяет достичь значительно более высокой точности результата, 
чем прямые методы идентификации, т. к. относительная ошибка оценки параметра немногим бо- 
лее 0,1 %. 
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РАКАМЕТКТС ТОЕМТТЕТСАТТОМ ОЕ САМ 5$У$ТЕМ МАТНЕМАТТСАЕ МООБЕ НАКОЕУ 
ОЕТЕВМТМАВЕЕ ТМУАВТАВЕЕ$ РОВ СА$ МАМ 5ЕСТТОМ$ 


В. А. Меудо!т, |. А. 2Пигамем, Е. \. Теегемема*, М. С. Егоо\у 
(Роп Зае Тесппка! УпмегКу, *УКа За е Тесйптса! Упмег$Ку) 


Тре аепийсавоп рго Мет оиЫШу о! зисй пагу а®еттта Ее рпуз/са! рагатеЁегх юг те дах тат зесйоп та- 
{Петайса! тодЕ! аз {те @5БиНоп ауегадеа а/опд те р/рете своп соейсепЕ апа {те гие теап епдЁТиауз 
ррете {етрегаиге ипеп изтд {те боептс тодЕе! мегэгоп [5 зи ед. Тре таетайса! аррагайи5 Гог 5исп 
АепЫйсавоп 15 деуеореа оп {те Баз/5 ОР 1е райегтеа ехрейтепЕ 

Кеуигога5: тат дах р/ре[те есйоп, САМ зу&ет, та{тетайса!/ тодЕ|/ $айопагу даз Йои,, рпузгса!/ соп5{апё 
рагатеис 'аепийсайоп. 
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